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D ' Abstract 

The structure function g\{x,Q 2 ) is analysed within the formalism based on 
unintegrated spin dependent parton distributions incorporating the LO Altarelli- 
Parisi evolution and the double ln 2 (l/x) resummation at low x. We quantita- 
tively examine possible role of the latter for the small x behaviour of the struc- 

"^j- ■ ture functions g\ of the nucleon within the region which may be relevant for the 

(N 

O 

structure function is dominated at low x by ladder diagrams the contribution 

0\ ' 

of the non-ladder bremsstrahlung terms is important for the singlet structure 

function. Predictions for the polarized gluon distribution AG(x,Q 2 ) at low x 



X 



O ,' are also given. 



1 Introduction 

Understanding of the small x behaviour of the spin dependent structure functions of 
the nucleon, where x is the Bjorken parameter is interesting both theoretically and 
phenomenologically. Present experimental measurements do not cover the very low 
values of x and so the knowledge of reliable extrapolation of the structure functions 
into this region is important for estimate of integrals which appear in the Bjorken and 
Ellis- Jaffe sum rules |J. Theoretical description of the structure function gf(x,Q 2 ) at 
low x is also extremely relevant for the possible polarised HERA measurements 0. 
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The purpose of the present paper is to explore the theoretical QCD expectations 
concerning the small x behaviour of the spin dependent structure functions taking into 
account the double ln 2 (l/x) resummation. The dominant contribution generating these 
double logarithmic terms is given by the ladder diagrams with the quark (antiquark) 
and gluon exchanges along the ladder. The very transparent way of resumming these 
terms is provided by the formalism of the unintegrated (spin dependent) parton dis- 
tributions which satisfy the corresponding integral equations. In this paper we extend 
this formalism so as to include the non-ladder bremsstrahlung terms. We argue that 
these terms can be easily included by adding the suitable higher order contributions 
to the kernels of the corresponding integral equations. We also incorporate within this 
scheme the complete LO Altarelli-Parisi evolution thus obtaining the unified system of 
equations which makes it possible to analyse simultaneously the parton distributions 
in the large and small x regions. In particular, this formalism allows us to extrapolate 
dynamically the spin dependent structure functions from the region of large and mod- 
erately small values of x, where they are constrained by the presently available data to 
the (very) small x domain which can possibly be probed at the polarized HERA. 

The content of our paper is as follows. In the next section we summarize the 
expectations of the Regge pole model for the small x behaviour of the spin dependent 
structure functions. We point out in particular potentially small magnitude (~ 0) of 
the corresponding Regge pole intercepts which control the small x behaviour of the 
spin dependent structure functions. We shall also remind that the Regge pole model 
expectations with the corresponding intercepts (~ 0) become unstable against the 
conventional QCD evolution which generates more singular behaviour at small x. In 
Section 3 we present the discussion of the double ln 2 (l/x) terms using the formalism of 
the unintegrated parton distributions. These distributions will be the basic quantities 
which will satisfy the corresponding integral equations generating the double ln 2 (l/x) 
resummation. We shall formulate these equations for the sum of ladder diagrams 
and extend the formalism by including the non-ladder bremsstrahlung terms. This 
will be done by a suitable modification of the kernel (s) of the corresponding integral 
equations. For fixed QCD coupling the equations will generate the complete small 
x behaviour of spin dependent structure functions. The formalism will be further 
extended by allowing the coupling constant to run and by including the complete LO 
Altarelli Parisi evolution. We shall discuss the analytic structure of the solution due 
to these modifications. In Sec. 4 we will present numerical solution of the integral 
equations starting from the simple semi-phenomenological parametrization of the non- 
perturbative part of the spin dependent parton distributions. Finally in Sec. 5 we give 
summary of our results. 



2 Regge pole model expectations for the small x 
behaviour of g\ and LO perturbative QCD effects 

The small x behaviour of spin dependent structure functions for fixed Q 2 reflects the 
high energy behaviour of the virtual Compton scattering (spin dependent) total cross- 
section with increasing total CM energy squared W 2 since W 2 = Q 2 (l/x — 1). This 
is, by definition, the Regge limit and so the Regge pole exchange picture |§ is there- 
fore quite appropriate for the theoretical description of this behaviour. Here as usual 
Q 2 = —q 2 where q is the four momentum transfer between the scattered leptons. The 
relevant Reggeons which describe the small x behaviour of the spin dependent structure 
functions are those which correspond to the axial vector mesons ||, |5j . 

The Regge pole model gives the following small x behaviour of the structure func- 
tions g[(x, Q 2 ) : 

g\(x J Q 2 )= ll (Q 2 )x-^°\ (1) 

where g\(x, Q 2 ) denote either singlet {gl{x, Q 2 ) = g{(x, Q 2 ) + g™(x, Q 2 )) or non-singlet 
(g™ s (x, Q 2 ) = <7i(x, Q 2 ) — g™(x, Q 2 )) combination of structure functions. 

In Eq. ([]]) «s,ns(0) denote the intercept of the Regge pole trajectory corresponding 
to the axial vector mesons with 7 = or / = 1 respectively. It is expected that 
«s,ns(0) < and that a s (0) ~ o> ns {$) i.e. the singlet spin dependent structure function 
is expected to have similar low x behaviour as the non-singlet one in the Regge pole 
model. 

Several of the Regge pole model expectations for spin dependent structure functions 
are modified by perturbative QCD effects. In particular the Regge behaviour ([]]) with 
«i(0) < becomes unstable against the QCD evolution which generates more singular 
behaviour than that given by Eq. ([[]) for ctj(O) < 0. In the LO approximation one gets: 



g^ix.Q 2 ) ~ exp[2JAP gg (0)aQ 2 )ln(l/x 



2f(*,Q 2 ) ~ exp[2y/x£Z(Q*)ln(l/x)], (2) 

where 



«*>-/;¥*£ (3) 



and 



_ AP qg (0) + AP gg (0) + ^/(AP M (0) - AP gg (0)) 2 + 4AP qg (0)AP gq (0) 

A - o I 4 ) 



with APjj(O) = APij(z = 0), where APy(z) denote the LO splitting functions describ- 
ing evolution of spin dependent parton densities. To be precise we have : 



AP qq (0) AP qg (0) 



N 2 -l 



ap(o) = 7 « - 7 « - = x. ; : , (5) 



AP OT (0) AP gg (0) 



N 2 - l 
N 




where N denotes number of colours, and Np denotes number of flavours. 



3 Double logarithmic ln 2 (l/x) corrections to g\ 

The LO (and NLO) QCD evolution which sums the leading (and next-to-leading) 
powers of ln(Q 2 /Qq) is incomplete at low x. In this region one should worry about 
another "large" logarithm which is ln(l/x) and resum its leading powers. In the spin 
independent case this is provided by the Balitzkij, Fadin, Kuraev, Lipatov (BFKL) 
equation [|7j which gives in the leading ln(l/x) approximation the following small x 
behaviour of the structure function Ff(x, Q 2 ) : 

F?(x,Q 2 ) ~ x - Xbfkl , (6) 

where the intercept of the BFKL pomeron Xbfkl is given in the leading order by the 
following formula: 

Xbfkl = 1 + —Un{2). (7) 

7T 

It has recently been pointed out that the spin dependent structure function gi at low x 
is dominated by the double logarithmic ln 2 (l/x) contributions i.e. by those terms of the 
perturbative expansion which correspond to the powers of ln 2 (l/x) at each order of the 
expansion [§, [| . Those contributions go beyond the LO or NLO order QCD evolution 
of polarized parton densities [[TIJ and in order to take them into account one has to 
include the resummed double ln 2 (l/x) terms in the coefficient and splitting functions 
rU [T2]. In the following we will discuss an alternative approach to the double ln 2 (l/x) 



resummation based on unintegrated distributions [13|, |TI 

To this aim we introduce the unintegrated (spin dependent) parton distributions 
fi(x',k 2 ) (i = u v ,d v ,u,d,s,g) where k 2 is the transverse momentum squared of the 
parton % and x' the longitudinal momentum fraction of the parent nucleon carried 
by a parton. Those functions will satisfy the corresponding linear integral equations 
generating the double ln 2 (l/x) resummation. 

The conventional (integrated) distributions Api(x,Q 2 ) are related in the following 



way to the unintegrated distributions fi(x', k 2 ): 

Ap t (x J Q 2 ) = Apf\x)+j k2 - w Mx' = x(l + —),k 2 ) } (8) 

where Ap\ (x) is the nonperturbative part of the distribution. The parameter k 2 , is 
the infrared cut-off which will be set equal to 1 GeV 2 . The nonperturbative part 
Ap\ (x) can be viewed upon as originating from the integral over non-perturbative 
region k 2 < k" 2 ,, i.e. 

Ap?\x)=f Q k2 °^-f t (x,k 2 ). (9) 

The spin dependent structure function g\ (x, Q 2 ) of the proton is related in a standard 
way to the (integrated) parton distributions: 



9 P i(x,Q 2 ) = 

(10) 



-(A«„(x, Q 2 ) + 2Au(x, Q 2 )) + -{Ad v {x, Q 2 ) + 2Au(x, Q 2 ) + 2As(x, Q 2 )) 

y y 



l 

2 

where Au v (x, Q 2 ) = Ap Uv (x, Q 2 ) etc. In what follows we assume that Aw = Ad and 
set the number of flavours Np to be equal to three. 

It is convenient to consider separately the valence quark distributions and the as- 
symetric part of the sea : 

fus(x',k 2 )=U(x',k 2 )-f s (x\k 2 ) (11) 

which do not couple with the gluons, and the singlet distribution : 

fc(x, k 2 ) = f Uv (x', k 2 ) + f dv (x', k 2 ) + Af n {x\ k 2 ) + 2f,(x', k 2 ) (12) 

together with the gluon distribution f g (x', k 2 ) which satisfy the coupled integral equa- 
tions. 

The full contribution to the double ln 2 (l/x) resummation comes from the ladder 
diagrams with quark and gluon exchanges along the ladder (see Figs. 1,2) and the 
non-ladder bremsstrahlung diagrams (Fig. 3). The latter ones are obtained from the 
ladder diagrams by adding to them soft bremsstrahlungs gluons or soft quarks [|[ || 
and they generate the infrared corrections to the ladder contribution. The sum of 
double logarithmic ln 2 (l/x) terms corresponding to ladder diagrams is generated by 
the following integral equations (see Figs.l, 2) : 



tv (MMM) i 



'"MMMtf 








Figure 1: Ladder diagram generating double logarithmic terms in the non-singlet spin 
structure function g± s . 



where j = u v , d v , us, and : 

Zn 



1 dz 



k I* dk' 2 
k' 2 



w>*)-w>*) + HtC%- 



AP qq (0)f s [j,k r 



AP 99 (0)/ s (^,fc' 2 



AP gs (0)/ 9 



AP gg (0)f 9 {-,k' 2 



(14) 



with APijiO) = APijiz = 0) given by Eq. (§). 

The variables k 2 (k' 2 ) denote the transverse momenta squared of the quarks (gluons) 
exchanged along the ladder, k$ is the infrared cut-off and the inhomogeneous terms 
f\ (a/, k 2 ) will be specified later. The integration limit k 2 / ' z follows from the require- 
ment that the virtuality of the quarks (gluons) exchanged along the ladder is controlled 
by the tranverse momentum squared. 

Equation (|13|) is similar to the corresponding equation in QED describing the double 
logarithmic resummation generated by ladder diagrams with fermion exchange |TB 



The problem of double logarithimc asymptotics in QCD in the non-singlet channels 



was also discussed in Ref. |T7|, [18], [19] 



Equations fll3|, |14|) generate singular power behaviour of the spin dependent parton 
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Figure 2: Ladder elements for a singlet case 



distributions and of the spin dependent structure functions g\ at small x i.e. 



' NS (x,Q 2 ) 



X 



-Ajvs 



^ ~x~ Xi 



9\ 

g*(x,Q 2 

AG{x,Q 2 )~ x - Xs , (15) 

where g± s — g\ — g™ and gf = g\-\-g"x respectively, and AG is the spin dependent gluon 
distribution. The behaviour of spin structure functions reflects the small x' behaviour of 
their unintegrated distributions. Exponents X^s s are given by the following formulas: 



A 



NS 



2tt 



AP M (0) 



(16) 



X, 



2tt 



Ao + 



(17) 



where Xq is given by Eq. (^). The power-like behaviour (|TJ) with the exponents Xns,s 
given by Eq. ( |ITD remains the leading small x behaviour of the structure functions pro- 
vided that their non-perturbative parts are less singular. This takes place if the latter 
are assumed to have the Regge pole like behaviour with the corresponding intercept (s) 
being near 0. The exponents X^s and Xs correspond to the leading singularities of the 
moment functions fi(u), k 2 ) ( Api(u>, Q 2 ) ) for % — NS(i = u v , d v , us) and i — E, i — g 
respectively in the uj plane. The moment functions ft(u, k 2 ) and Api(u, Q 2 ) are defined 




c# 




a) 



b) 




Figure 3: Non-ladder contributions containing a), b) soft gluon or c) soft quark attached 
to the quark-quark scattering amplitude 

in a conventional way : 



fi(u,k')= / dx'x'"- l fi{x\k 2 ), 
Jo 



Ap i (u,Q 2 )= f 1 dxx^Ap^Q 2 ). 

Jo 

Ladder equations flIB|), ( [HP in u— space take then the form : 



Uw,k 2 ) = f( w ,^) + ^AP M (0) 



«„*^(- 4 ' 2 » + ri?(^) *(»•*") 



f s ( U ,* ! ) = 4°'(a,,t 2 ) + *M 



* ^(«.* a )+jT^(^) «.(-.*-) 



where i = it„, <L, us, and a s is defined as : 



O-yj u*-y j U/O, 



«.s 



and 



frfw,*: 2 } 



M = AP(0) 



_ «5 

= 2tt' 

h(oj,k 2 ) 
f 9 (",k 2 ) 

AP qq (0) AP qg (0) 
AP gq (0) AP gg (Q) 



,(19) 
(20) 

(21) 



(22) 






Figure 4: Infrared evolution equation for partial wave matrix F$ 

The inhomogeneous terms fl and fg are related to the moment functions Ap\ (a>) 
of the nonperturbative parts Ap\ (x) and Ap^ (x) of the integrated distributions in 
the following way : 



m 



fr {<*,*?) = a 



* AP w (0) 



*(°) 



w 



AiTM, 



ff(u,,fc 2 ) 



_ Mo A _(o)/ x 



UJ 



For fixed coupling a s one may obtain analytic solution of (|i9|), ( gO 



fi(u,k 2 



Ri(a s ,uj) 



'k_ 



2\ 7ats("s,w) 



f5(-,^)=(| 



2\ 7s( Q =' w ) 



R s (a a ,u;) 



(23) 



(24) 



(25) 



where Jns, 7 s represent anomalous dimensions for the non-singlet and singlet case 
respectively, and are equal to : 



UJ - 


-^ 


-4AP qq (0)a s 






2 


UJ - 


-Vu 2 


- 4M a s 



j NS {uJ,a s ) -- 

and Ri(a s ,uj) and Rs(a s ,w) are given by the following equations : 

Ri(a s ,uj) = >y N s(a s ,uj)Apf\uj) 
~R s (a s ,uj) = -y s (a SJ uj)Ap ( g\uj). 



(26) 
(27) 



(2* 
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It should be reminded that for the singlet case only the eigenvalues of "y s contribute 
to the final solution. The solution takes then the explicit form : 



'£,2\ 7 + (^.as) fk 2 ^ 1 ^'"^ 



h, g (u, k 2 ) = R+ ig (u, a s , ) [ ^ I + i^>, a,) 1^1 , (29) 

where Xq, Xq denote eigenvalues of of 7 5 : 

7 ± (w,a,) = * — , (30) 

and \q (the eigenvalue of M ) is defined by equation (f|) while A^~ reads: 



_ AP M (0) + AP gg (0) - ^(AP qq (0) - AP gg (0)) 2 + AAP qg (0)AP gq (0) 
A o g ' ' ' 

It can be shown that the moment functions Api(u),Q 2 ) and Ap s (u,Q 2 ) of the spin 
dependent distributions have the familiar RG form (for the fixed coupling) : 

//02\ 7ivs(<5s,w) 

Ap t (Lu,Q 2 ) =Ri(a a ,u) f ^-J +0(k 2 /Q 2 ) 



Ap,(^.Q-) ■ (tt! " Rx-ln^ + OU-oVQ-). (32) 



where 



r(w) 

ts /_ v r(7 g (a s , a;) + l)r(o; - -f s (a s , top) (0 ) , , , . 

R 5 (a s ,^) = — Apy(w). (33) 



To derive fl32 ) we notice that equation (|8|) implies the following relation between the 



moment functions Api(u,Q 2 ) and fi(u,k 2 ): 

Aft( w> Q 2 ) = Apf ( W ) + jT ^ f l + |V /<(", fe 2 )- (34) 



Substituting expressions (]24]) and ( p5|) into ([34]) and taking into account equation (p8|) 
we obtain equation fl3~2"|). 

Equations fl26|), (|27|) imply that the anomalous dimensions 7^5(0;, a s ) and 7 ± (w, a s ) 



have the branch point singularities in to plane located at u — 2Ja s AP qq (0) and at 



uj = 2Ja s X Q respectively. The small x behaviour of the parton distributions and of 
the structure functions which is given by equations (|HJ) with the exponents X^s an d Xs 
defined by equations (pT6|) and (|T7| ) just reflects the fact that this behaviour is controlled 



by the leading singularities in the u plane. 
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It may be seen from equations (^,^) that the branch point singularities of anoma- 
lous dimensions are also present in the prefactors Ri(a s ,u>) and Rg(a s ,w) and so the 
low x behaviour of spin dependent structure functions given by equations fll5|) is ex- 
pected to hold for arbitrary values of the scale Q 2 . 

Complete sum of double logarithmic ln 2 {l/x) terms does also include the non- 
ladder bremsstrahlung contributions besides the ladder ones. Expressions (|IED , (|TTD 
corresponding to the ladder diagrams are therefore only approximate. The method of 
implementing the non-ladder bremsstrahlung corrections into the double logarithmic 



resummation was originally proposed by Kirschner, Lipatov [ I7| , and applied for the 



case of DIS small x asymptotics in || [|, ]TTJ. Below we show how to implement 
these terms within our formalism. In order to include the non-ladder corrections to 
the double logarithmic asymptotics we use the infrared evolution equations derived in 



17| . The evolution equations for singlet partial waves F , F$ (see Fig. 4 ) : 

(35) 



0,8 



ei(0,8) f(0,8) 

m 19 

F(0,8) p(0,8) 
99 99 



have the form 



7T a s , , 4a, 



F (u>,a s ) = ^^M - ^F 8 (u,a s )G + -^F 2 (uj,a s ) (36) 

u or 8ti z uj 

F 8 a;,a a = -M 8 + ^—F 8 (u,a s ) + -^F 2 8 (u,a s ), 37 

u uj duj 8tt z uj 

where M (cf. ^), M 8 are splitting functions matrices in colour singlet and octet 

t— channel, and Go contains colour factors resulting from attaching the soft gluon to 

external legs of the scattering amplitude : 

M s = 




Go = 2 : : T ■ (a 



Non-singlet contribution comes from Fq q partial wave, as expected. Only the partial 
waves Fo are contributing to the anomalous dimensions. An explicit relation between 
partial waves and anomalous dimension reads : 

Fo = 8vr 2 7 f^ (39) 

Fj$ = 8^§f (40) 

and it depends already on resummed anomalous dimensions ^ RES with all infrared non- 
ladder corrections included. Equation (|36|) for F may be solved analytically generating 
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the following expressions for anomalous dimensions : 



7#s («.><*>) 



7s ES (a s ,uj) 



u> 



U) 



1 Aaj AP qq (0) (F 8 MG 



'gy 



2vr 2 ^ 2 



\ 



4a s /M F 8 (u)Gq 
2ir 2 u 2 



oj \ UJ 



(41) 
(42) 



Comparing the expressions for anomalous dimensions with and without the non-ladder 
contribution (|26|), fl2~T|), (f41"D, ( |42| ) one may notice that the bremsstrahlung contribution 
adds an additional term under the square root of 7's : 



4a, Mr 



UJ UJ 



\ 



4a, /M n F 8 MG r 



UJ \ UJ 



2ix 2 uj 



2,. ,2 



(43) 



The same anomalous dimensions would be obtained if we modified the kernels of equa- 
tions (|]JD and (H) by setting a s (^ - F ^g° ) in place of a s ^. The modified 
equations then read: 



_ f AP(0) FsHGq 

f s ( W ,fc 2 ) = f£V,fc 2 ) 

Mq F 8 (^)G 

W 27T 2 W 2 



'/'/ 



dk'' 



l>2 fi(u,k' 2 ) + 

k 2 K Jk 2 



00 WZ-'2 / Z.2 \ u 



dk 

"fc 72 Vfe 



/i(w,fc' 2 ) 



dfc' : 



00 WJ,'2 



dfc' 2 / k 2 \ 



;j ^(^' 2 )+y fe2 ^(^J ^^ 



(44) 



(45) 



The corresponding equations for the functions fi(x', k 2 ) and for for f^ have the form : 

- 1 dz r k2 / z dk' 2 „ fx' 



f J (x',k 2 ) = ff( X ',k 2 ) + a s AP gg (0)/jf/J /Z ^(^^) 

.„ \ r k2 dk' 2 . fx' 



O.s 



O.s 



dz 



1 z 



UJ^ 



(*) 





Gc 

2vr 2 



-/a lk ° 



k' 2 



*(t^) 



1 ^ /- fc2 /* dfc' 2 
"fc 72 " 



a;' 2; Jfc2 



W* 



A; 2 * ) 2vr 2 



I /i(pfe /2 ),(46) 

99 



where j = u v , d v , us, and : 

Hx\k 2 ) = if\x\k 2 ) + a s M f l -f 

Jx' Z Jk 



1 dz r k2 / z dk' 2 



f.<? -,* 



/2 



n, s 



a. 



dz 



' z 



a; 



feg fc' 2 ° V z 

G /" fc2 <iA: /2 r ( x ,,, 



m^L^Hp* 



1 dz r k2 / z dk' 2 



z Jk 2 



k' 2 



UJ^ 



C^u( x -y> 



k 2 2tt 2 



(47) 



12 



where 



Ea 



denotes the inverse Mellin transform of ^4 



d + io ° du _ w F 8 {u) 



u^ 



5— ioo Z7TZ 



<jJ z 



(4£ 



with the integration contour located to the right of the singularities of the function 

F 8 H 



To obtain F 8 we have to solve Eq. flT7|). The exact solution presented in R [| may 
be expressed in terms of parabolic cylinder functions : 



/I 



du 



(49) 



where / 8 are two eigenvalues of matrix Fg(u;) determined in a basis of eigenvectors of 

M 8 , and : 

A^ 



z = — , u 
u 



Na s , p± = — , 



(50) 



where A^ denote eigenvalues of matrix M 8 . 

For a general solution fl4*9| ) the inverse Mellin transform of 2 does not exist in 
the analytic form. However, it was checked in || that approximate form of F$(uj) 
determined in the large N limit is a good approximation for fixed a s . The large N 
expansion proposed in H reduces Fg to the Born term. In our case it implies : 



iBorn 



(to) ~ 8n 2 a, 



M 8 



u 



The inverse Melin transform then reads 

~T?Born 



U* 



z) = 47Ta,M s Zn r(z) 



(51) 



(52) 



The different large N approach which we propose here uses the approximate form of 
M , M 8 , G (H), © : 



Mr 



N/2 

N AN 



M 8 





iV 2N 



1 G 



N/2 
N 



(53) 



where we have neglected all terms of the order less than O(N). Then the two compo- 
nents of F 8 simplify : 



ft " ' ' (-') 



fs 



,large N 
,large N . 



IQtt 2 



UUn 



u* — u>, 



2> 



u) 



(54) 
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/■+, large N 

and the inverse Mellin transform of ^ — 5 — reads : 



n 



-,largeN 



uJ 2 



ir^z" + z'^ - 2). (55) 



We have checked that both approaches give similar evolution of polarized structure 
function gi and gluon distribution. Moreover, the branch point singularities dominating 
the small x behaviour of g± s and gf behave also in similar way. For the non-singlet 
case the leading singularity may be determined from the Eq. (pfl| ) as : 



45, AP qq (0) _ (F 8 (u)G ) qq 
u> \ uj 2tc 2 lu 2 



1 _~,, -~ wv -, _ v- y "»/ m = Q (56) 

and for the ladder case it reads (cf. fllfiD ) 



X NS = 2^pS. « 0.39, (57) 

where it was assumed that a s = 0.18, N = 3, Np = 3. Including the non-ladder 
resummation gives : 



A „, flm „ J^ 1 + ^W 



A 7VS,/ar 9 eV _ 2d —a s ^ 0A2. (5* 

The singlet dominating singularity fulfills the relation (cf. (f42|)) : 

4a s /M F 8 HG ' 



d e t (l-^£(^-ii^ =0 (59) 



and for the ladder case it may be easily determined as : 

X s = 2\/\£a s & 1.13 (60) 

(a s = 0.18, N = 3, Np = 3). The non-ladder resummation changes the singularity 
point into : 

X S,Born _ LQ1 
X S,lar 9 eN _ j_ q ? (61) 

One may notice that for both singlet and non-singlet case the influence of the non- 
ladder resummation on the singularities is of the order 10%. Since the differences 
between Born and large N approaches are very small, in the following discussion we 
restrict to the Born approximation for Fg. 
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4 Unified evolution equation 

In order to make the quantitative predictions one has to constrain the structure func- 
tions by the existing data at large and moderately small values of x. For such values of 
x however the equations ( pEop and (17) are inaccurate. In this region one should use the 
conventional Altarelli-Parisi equations with complete splitting functions APij(z) and 
not restrict oneself to the effect generated only by their z — > part. Following refs. 
15| we do therefore extend equations (|46l , |4"7|) and add to their right hand side the 



contributions coming from the remaining parts of the splitting functions APij(z). We 
also allow the coupling a s to run setting k 2 as the relevant scale. In this way we obtain 
unified system of equations which contain both the complete LO Altarelli-Parisi evo- 
lution and the double logarithmic ln 2 (l/x) effects at low x. The corresponding system 
of equations reads : 



fk(x',k 2 ) = fl U) (x',k 2 ) + 



(0), 



+ 



+ 



a s {k 2 )4 f 1 dz r k2 i 
2n 3 Jx' z Jkl 


z dk' 2 fx' 
k' 2 h \z ' 


k' 2 \ 




( Ladder ) 




a s (k 2 ) r k2 dk! 2 4 r 
2tt Jkl k' 2 3Jx 


idz(z + z 2 )f k (^,k' 2 )- 


2zf k (x',k' 2 ) 


' z \ — z 




a s (k 2 ) r k2 dk' 2 
2tt Jkl k' 2 


2+-ln(l-x') 


fk(x',k' 2 ) 




( A-P ) 




a s (k 2 ) j^dzf' 
2ir J x ' z y 


F 8 " 
w 2 


G \ f dk' 2 f 
{z) 2n 2 Jkl k' 2 ^{ 

/ qq 


7-* B ) 


a s (k 2 ) r l dz r k ' 2 / z dk' 2 ( 
2tt Jx> z Jk? k' 2 y 


F 8 " 

UJ 2 


fk' 2 \G N 
\k 2Z j2n 2 / 


QQ 


(Non 


-la 


dder) 











(62) 



where k = u v , d v , us, 

f f > &\ f(°)( ' i, 2 \ , a "^ f ldz f k2/zdk ' 2 4 f fx' ,, 2 
Mx,k) = f^( X ,k) + -j-J^-J -^^[j.k 



+ 



+ 



+ 



a s (k 2 ) f 1 dz f k ' 2 / z Ar dk'' 



2tt 



V7L' "'yrSA--:,* 



( Ladder ) 

n s (k 2 ) [# dk 12 4 /"I dz (Z + Z 2 )M^, k' 2 ) - 2zfe(x', k' 2 ) 

\-z 

Mx',k' 2 ) 



2tt J^ k' 2 3 J X ' z 



a s (k 2 ) 



a. 



2tt 



dk 



.12 



2tt Jkl k' 2 



2 + -ln(l-x') 



^ dk' 2 f l dz 



f dk' z f 1 dz n AT , ,x' ,«, 
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( A-P 



a s (k 2 ) f l dz 
2tt J x ' z 



2tt 



(*) 



G \ f dk'\ (x> 

/ qq 



•2 



a s {k 2 ) f 1 dz f k l z dk 



12 



z Jk 2 k' 2 \ 
(Non — ladder) 



u' 







k 2Z )2ir 2 



/', 



"9 



(63) 



-,k' 2 



f g (x',k 2 )=6 \x',k 2 ) + 



+ 



+ 



+ 



+ 



a s {k 2 ) f 1 dz_ r k2 / z dk^8 fx[_ 2 
2tt L z hi k' 2 3 /E U' 



n.s 



(fc 2 ) f 1 dz f k2 / z dk 



2n 

a s {k 2 \ 
2vr 

a s {k 2 ) 
2vr 

a s {k 2 ) 
2vr 

a s (k 2 ) 
2ir 



i z k 2 



/2 



k' 2 



Wo 



A: 



12 



( Ladder ) 
f k ' 2 dk' 2 f l dz . 4, „ fx' , 



k2 dk' 2 rl 



k 2 k' 2 J x > : 



dz 
— oz 



f 9 (^,k' 2 )-f g (x',k' 2 ) 



VA-^ 12 



k2 dk' 2 



2 k' 2 





( 

1 dz 

■' z 



11 N F 



+ 6ln(l-x') 



' lJ2\ 



Jg\ x >k 



A-P 

F 8 



) 



(*) 



gA 



k2 dk' 2 



2fT 2 ) J k 2 k' 2h \z' k 
/ 00 o 



.00 



/2 



(64) 



a s {k 2 ) r 1 dz r k2 / z dk' 2 ( 

2vr L T Jk 2 ~k 1T [ 

(Non — ladder) 



UJ^ 



VMGc/ 

k 2 Z ) 2vr 2 



h 



Ml 



J2 



In equations (|62|), (p3|), (|6J) we group separately terms corresponding to ladder dia- 
gram contributions to the double ln 2 (l/x) resummation, contributions from the non- 
singular parts of the Altarelli-Parisi splitting functions and finally contributions from 
the non-ladder bremmstrahlung diagrams. We label those three contributions as "lad- 
der", "AP" and "non- ladder" respectively. 

The inhomogeneous terms fl (x',k 2 ) are expressed in terms of the input (inte- 
grated) parton distributions and are the same as in the case of the LO Altarelli Parisi 
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/JV,* 2 ) 



a s {k 2 ) 4 /-i dz (1 + z 2 )Ap k 0) {^) - 2zA Pk °\x' 



,(o), 



2it 3 Jx' z 
a s {k 2 



2ti 



2 + -ln{\ - x 1 ) 



1-z 

Ap k °\x') 



(65) 
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I /C — ^V) d"Vl US ) ^ 



,(0) 



f^{x\k 2 ) 



a s (k 



4 /-i dz (1 + ^ 2 )Apg ] (|) - 2zApx 

2rr 3 



(°), 



.?: 
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1-2 



+ ^(2 + |ln(l-*0)Ap< 3 V) 



27T Jx' Z 



2z)Apf{ 



fi°\x',k 2 



a s (k 2 ) 4 yi cfe 



27T 3 A' 2 



-(2-*)Ap^ 



a s (A; 2 ) ,11 JV F „, .„ 



Co), 



;r 



(66) 



Q! g (fe 2 ) /^dz 

2-7T 7a;' 2 



^bW +[1 . 2;|Af ,„ ( < 

1 — y z 



Equations (0), (|63]) , fl64] ) together with fl65[), (|66|) and (|[) reduce to the LO Altarelli- 
Parisi evolution equations with the starting (integrated) distributions Ap®(x) after we 
set the upper integration limit over dk' 2 equal to k 2 in all terms in equations (|62|), (j63|) , 
(|64]), neglect the higher order terms in the kernels, and set Q 2 in place of W 2 as the 
upper integration limit of the integral in Eq. (^). 

The presence of the running coupling changes the singularity structure of the so- 
lution turning the branch point singularities into the infinite number of poles whose 
position depends upon the magnitude of the cut-off k 2 |19| . Apparent branch point 
singularities are present if we adopt the semiclassical approximation to the solution of 



equations (|62"D, (|63"D, (|64"D with the running coupling. In this approximation we just 
recover the RG structure with the running coupling i.e. 

„2\ 7ivs(^,as) 






exp 



k 2 dk 



12 



kl k 



/2 



■~f NS (u},a s (k )) 



(67) 



and similarly for the singlet part. 

Introducing the Altarelli-Parisi kernel into the double logaritmic evolution equations 
(^).(^) changes the anomalous dimensions 7 (|4"T|), (f£2|). They take then the form : 

u + H qq (uo) - J(uo - H qq (u)) 2 - AP qq (u)uj 



-f NS (uJ,a s ) 



IS (W. a s 



u + H(w) - J(u - B.(lo)) 2 - 4~P(lu)uj 



(6S 



where H(u) denotes the moment representation of non-singular part of the Altarelli- 
Parisi kernel, and P(oS) denotes the moment of double logaritmic kernel i. e. for ladder 
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case P{oj) = as ■ We investigated the influence of the Altarelli-Parisi kernel on the 
behaviour of leading sigularities for the non-singlet and gluonic sector (see Tab. 1), 
assuming for simplicity : 

HM 



PM 



H(w = 0) 
Pi P 2 



uj 



U c 



(69) 



where Pi(w) = a s M and P 2 (w) = 4<^M 8 G . 



In Table 1 we summarize numerical results concerning the position of leading branch 
points for the non-singlet and singlet case obtained in different approximations of the 
kernel i.e. for the pure double logarithmic approximation generated by ladder diagrams 
alone (first column), for the complete double logarithmic approximation (second col- 
umn) and finally for the complete double logaritmic approximation supplemented by 
the non-singular part(s) od the Altarelli-Parisi splitting functions (third and fourth 
column). One may notice that the position of the branch point for the non-singlet 
part is to a very good accuracy determined by the double logarithmic approximation 
generated by ladder diagrams alone, and the non-ladder contributions and the non- 
singular parts of the AP kernel do not play important role. The latter terms are 
however important in the singlet case, and in particular the non-singular part of the 
Altarelli-Parisi splitting function significantly reduce the magnitude of the position of 
the branch point singularity. This means that in the singlet case the corrections to 
the double logarithmic approximation may be expected to be very important. This 
effect will be quantified in the next Section where we present results of the numerical 
solution(s) of equations (|62|) , d63|) , (|64|). 

Table 1. Leading singularities for the non-singlet and gluonic part of g\ determined for 
a s = 0.18, N = 3,n f = 3. 





ladder 
H = 0,P 2 = 


+non — ladder 
H = 




ladder + (A - P) 
P 2 = 


+non — ladder 
+(A-P) 


X NS 


0.39 


0.41 




0.41 


0.43 


X s 


1.17 


1.08 




0.96 


0.78 



4.1 Numerical results 

We solved equations (|62|), (|63|), (|64]) assuming the following simple parametrisation of 
the input distributions: 
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Ap, 



C°), 



X 



NAl-x)^ 



(70) 



with f] Uv = rj^ = 3, r] a = t] s = 7 and rj g = 5. The normalisation constants AT, were 
determined by imposing the Bjorken sum-rule for Au^ — Ag^ and requiring that 
the first moments of all other distributions are the same as those determined from the 
recent QCD analysis [pl| . All distributions Ap\ (x) behave as x° in the limit x — > 
that corresponds to the implicit assumption that the Regge poles which should control 
the small x behaviour of g[ have their intercept equal to 0. 
10 
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Figure 5: Non-singlet part of the proton spin structure function gi(x, Q 2 ) as a function 
of x for Q 2 =10 GeV 2 . Solid line corresponds to the calculations which contain the full 
ln 2 (l/x) resummation with both bremsstrahlung corrections and Altarelli-Parisi kernel 
included, dashed line represents the ladder ln 2 (l/x) resummation with Altarelli-Parisi 
kernel included, a dotted line shows the pure Altarelli-Parisi evolution, and a thin solid 
one describes the input non-perturbative part g[ . 



It was checked that the parametrisation (ffD|) combined with equations (H), (|l0l),(|62l), 
gives reasonable description of the recent SMC data on g^ s (x,Q 2 ) and on 
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Figure 6: The structure function c/i(x, Q 2 ) for Q 2 = lOGeV 2 plotted as the function of 
x. Solid line corresponds to the calculations which contain the full ln 2 (l/x) resumma- 
tion with both bremsstrahlung corrections and Altarelli-Parisi kernel included, dashed 
line represents the ladder ln 2 (l/x) resummation with Altarelli-Parisi kernel included, 
a dotted line shows the pure Altarelli-Parisi evolution, and a thin solid one describes 
the input non-perturbative part g^ . 



gi(x,Q 2 ) f22| . In Fig. 5 we present the nonsinglet part of gi(x, Q 2 ) for Q 2 = lOGeV 2 

and confront 



in the small x region jl3| . We show predictions based on equations (| 
them with the results obtained from the solution of the LO Altarelli-Parisi evolution 
equations with the input distributions at Ql = lGeV 2 given by equations (|70|) . We 
also show results based on equations similar to flG2| ) in which we have removed the 
bremsstrahlung contributions to the kernel. We may see from this figure that the double 
logarithmic contributions are very important at low x, and that they are reasonably well 
described by the contribution of ladder diagrams. We also plot the nonperturbative 



part of the non-singlet distribution g x 
vector coupling. 



NS(0) 



X 



W6(l 



x 



where gA is the axial 



In Fig. 6 we again confront predictions for g\ (x, Q 2 ) at Q 2 = lOGeV 2 based on 
equations (|| |63|, 0) with those based on the LO Altarelli-Parisi evolution equations, 
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Figure 7: The spin dependent gluon distribution AG(x, Q 2 ) for Q 2 = lOGeV 2 plotted 
as the function of x. Solid line corresponds to the calculations which contain the full 
ln 2 (l/x) resummation with both bremsstrahlung corrections and Altarelli-Parisi kernel 
included, dashed line represents the ladder ln 2 (l/x) resummation with Altarelli-Parisi 
kernel included, a dotted line shows the pure Altarelli-Parisi evolution, and a thin solid 
one describes the input non-perturbative part AG' '. 

and with the results based on equations similar to (|63| , |64[) in which we have removed 
the bremsstrahlung contributions to the kernel. In the region of very low values of 
x the dominant contribution comes from the singlet component of g\. We see from 
this figure that the contribution of the bremsstrahlung term is very important and 
significantly slows down the increase of g\ with decreasing x. The structure function 
gi(x, Q 2 ) which contains effects of the double ln 2 {l/x) resummation begins to differ 
from that calculated within the LO Altarelli Parisi equations already for x ~ 10~ 3 . 

In Fig. 7 we show the spin dependent gluon distribution which contains effects of 
the double ln 2 (l/x) resummation and confront it with that which was obtained from 
the LO Altarelli-Parisi equations. It can be seen that the ladder resummation exhibits 
characteristic x~ Xs behaviour with Xg ~ 1. Similar behaviour is also exhibited by the 
structure function g\ (x, Q 2 ) itself. The contribution of the bremsstrahlung term is very 
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important and significantly slows down the increase of AG with decreasing x. 



5 Conclusions 

To sum up we have presented theoretical expectations for the low x behaviour of the 
spin dependent structure function gi(x, Q 2 ) which follows from the resummation of the 
double ln 2 {l/x) terms. We have also presented results of the analysis of the "unified" 
equations which contain the LO Altarelli Parisi evolution and the double ln 2 (l/x) 
effects at low x. We based our calculation on a formalism of the unintegrated spin 
dependent distributions which satisfied the corresponding integral equations. This for- 
malism made it possible to make an insight into physical origin of the double ln 2 (l/x) 
resummation. The structure of the corresponding equations is similar to the conven- 
tional LO Altarelli - Parisi evolution equations with extended kernels to account for 
the bremsstrahlung contributions. Very important difference however is the absence 
of transverse momenta ordering along the chain (let us remind that the LO Altarelli- 
Parisi evolution corresponds to ladder diagrams with ordered transverse momenta). 
This should have important implications for the structure of the final state in polar- 
ized deep inelastic scattering at low x, similarity to the case of unpolarized DIS |2B| . 
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